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Steady state entanglement of two coupled qubits 
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The maximum entanglement between two coupled qubits in the steady state under two independent incoherent 
sources of excitation is reported. Asymmetric configurations where one qubit is excited while the other one 
dissipates the excitation are optimal for entanglement, reaching values three times larger than with thermal 
sources. The reason is the purification of the steady state mixture (that includes a Bell state) thanks to the 
saturation of the pumped qubit. Photon antibunching between the cross emission of the qubits can be used to 
experimentally evidence the large degrees of entanglement. 
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I. INTRODUCTION 

A global state of a composite system is entangled when it 
cannot be written as a product of the states of the individual 
systems ^ . This is the basic quantum-mechanical property, 
with no classical analog, for quantum information technolo- 
gies 12, S]. Two coupled qubits, or two-level systems with 
ground | g) and excited | e) states, is the smallest and simplest 
composite system that can display entanglement. It is, there- 
fore, the most suitable model to investigate its creation and 
processing as well as how environmental noise and decoher- 
ence brought by spontaneous decay and the external excitation 
affects it |4|], which is a key point for quantum applications. 

Two qubits can form four independent maximally entan- 
gled states, the so-called Bell states: \(^±) = (\gg)± \ee))/\/2 
and |v/±) = {\eg) ^ \ge))/\/2. The last two are affected by 
a possible coupling between the qubits, and are also known 
in the atomic literature as the symmetric and antisymmetric 
collective states ^ . The formation and degradation of such 
states when subjected to spontaneous emission has been the 
object of much recent research |[a-[sl], focusing on the preser- 
vation of entanglement into decoherence-free subspaces and 
taking advantage of the collective damping or effective cou- 
pling created between the qubits by interaction with common 
reservoirs |0-[2fl]- 

The idea of environmentally induced entanglement has also 
been applied to the case where the two qubit interaction is me- 
diated by a cavity mode (harmonic oscillator) which is excited 
by white noise (a thermal reservoir) [211 , borrowing the idea 
from Ref. ll22n where, on the contrary, entanglement is en- 
hanced between two harmonic modes by mediation of a two- 
level system excited by white noise. In both cases, entangle- 
ment may survive in a steady state that is not the vacuum, but 
is very small (< 0.4%). A two-level system has also been 
proposed as a mediator (or coupler), to build entanglement 
between qubits ll23[ [24ll. 

Another possibility, close to the one addressed in the 
present text, is to consider two qubits already coupled, whose 
entanglement builds in the steady state despite dissipation 
and decoherence from two independent environments 11251 - 
[33n . Entanglement, being essentially a property that requires 
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great purity of the state, is very sensitive to such decoher- 
ence. It is therefore important to look for optimization. Hav- 
ing a high degree of entanglement in the steady state means 
that it is robust and independent of the initial state, it remains 
stored forever in our system like a quantum battery of entan- 
glement 132]. 

The coupling between the qubits can have different physical 
origins depending on the realization 13J, e.g., Rydberg atoms 
couple through dipole-dipole interaction, weaker in the case 
of cold atoms p?], and so do excitons in single quantum dots 
or molecules |35-37]; superconducting qubits couple through 
mutual inductance Il38ll . Moreover, in all these implementa- 
tions, the coupling can appear effectively through the virtual 
mediation of a coupler (a cavity or a wire mode) in the disper- 
sive limit, in which case it is given by geff ^ ^c/^c, where 
Gq is the coupling of the qubits to the coupler and Ac the 
energy detuning to the qubits (considered much larger than 
the coupling) Il39l - [4in . This scheme requires, for instance, 
placing the qubits into a cavity where the cavity mode acts as 
the coupler. One can take advantage of the QED techniques 
while obtaining an effective coupling essentially insensitive to 
the cavity decay and thermal fluctuations. The effective cou- 
pling between two Rydberg atoms through virtual photon ex- 
change, while crossing a nonresonant cavity, was achieved in 
2001 142]. The final entangled state could be controlled by ad- 
justing the atom-cavity detuning. A similar effective coupling 
was obtained between two superconducting qubits on opposite 
sides of a chip using microwave photons confined in a trans- 
mission line cavity 11431] . The cavity was also used to perform 
multiplexed control and measurement of both qubit states. Ef- 
fective coupling between two distant quantum dots embedded 
in a microcavity has also been recently achieved 14^ SSfl . 

Taking for granted that the two qubits are coupled, we 
center our attention on the situation where the qubits are 
also in contact with two independent excitation sources. Xu 
and Li 11251] found that with two equally intense white-noise 
sources at the same temperature, no entanglement can be 
formed in the steady state. However, if only one qubit 
was subjected to a finite temperature source, some entan- 
glement could be achieved. They also pointed out that the 
steady state entanglement exhibits a typical double stochastic- 
resonance as a function of the decoherence parameters of both 
qubits 1 22] . They found better but still small degrees of en- 
tanglement (< 4%) and did not deepen on its origin, but their 
results show that an asymmetric flow of excitation through the 
qubits is beneficial for entanglement. Other authors, who did 
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FIG. 1 : Scheme of the two coupled qubits or two-level systems (a) 
and their energy levels (b), with coupling (g), pumping (Pi) and decay 
parameters (yi). 



not consider unequal sources of excitation (having to recur to 
other mechanisms for entanglement generation), explored, on 
the other hand, configurations that are out of thermal equilib- 
rium, where the excitation of the qubits is not necessarily pro- 
duced by thermal sources but by more general p rocesses able 
of inverting the qubits population 1271 IZSl ^M- Two qubits 
may undergo dissipation and pure dephasing but also an ex- 
ternally controllable and independent (in general) continuous 
pumping that can have great impact on the strong coupling 
reached in the steady state ll46l l47ll. 

In the present text, I put together different elements that 
have been addressed separately in previous studies of envi- 
ronmentally induced entanglement: direct coupling between 
the qubits and independent and different kinds of reservoirs 
that are not necessarily of a thermal nature. I give a com- 
plete picture of entanglement and its origin in the steady state 
of such a general system. As a result, I find a configuration 
where entanglement is significantly enhanced (31%), that is, 
more than three times as compared to the best thermal case 
and with a much better purity. I show that it can be evidenced 
by the antibunching of the two qubits cross emission. 

The rest of this paper is organized as follows. In Sec. [Ill I 
introduce a theoretical model to describe two coupled qubits 
with decay, incoherent pumping and pure dephasing, and a 
quantity to quantify the degree of entanglement between them, 
the concurrence. In Sec.|llll I discuss different entangled con- 
figurations and optimize the concurrence for the most suitable 
one: one qubit is excited while the other dissipates the ex- 
citation. This is compared with the thermal counterpart. In 
Sec. HVl I show how a strong antibunching between the two 
qubits emissions, is linked with high degrees of entanglement 
and I propose this effect as an indication of entanglement. In 
Sec.|Vl I study the effect of pure dephasing on entanglement. 
In Sec.|Vll I present the conclusions. 



II. THEORETICAL MODEL 

Let us consider two qubits or two-level systems (/ = 1,2), 
with lowering operators O";, frequencies cOf and coupled with 
strength g. Without loss of generality, we take the energy of 
the first qubit as a reference (coi =0), from which the other 



one is detuned by a small quantity A= coi 
sponding Hamiltonian reads: 



CO2. The corre- 



H = -AgIg2 + g(cr/cr2 + crj'cTi) , 



(1) 



Each qubit can be in the ground (\g)) or excited (|^)) states, 
whose direct product produces a Hilbert space of dimension 4 
(see Fig.ID: {|0) = \gg), |1> = \eg), |2) = \ge), |3) = \ee)}. 
The qubits are in contact with different kinds of environments, 
that provide or dissipate excitation (at rates Pi and Yi^ respec- 
tively) in an incoherent continuous way. Other interactions 
may purely bring dephasing to the coherent dynamics (at rates 
Yf). These processes eventually drive any pure state into a 
statistical mixture of all possible states. A density matrix, p, 
properly describes the evolution of such a system. The gen- 
eral master equation we consider has the standard Liouvillian 
formS: 



a,p = /[p,if] + £ 



i=i 



Yi 



°2-(T, 



2 <^l 



yf 



^ 



(2) 



with the corresponding Lindblad terms for the incoherent pro- 
cesses {^op = 20pO^ - O^Op - pO^O). If the two qubits 
shared a common environment, the Lindblad terms would 
share a single expression ^j in terms of the collective op- 
erator: / = CTi + (72. Such collective terms are sources of 
entanglement, as explained in the introduction. In this text, 
I investigate the steady state of a system where they are not 
present, solving exactly the equation dtp =0. 

In order to spell out the nature of the reservoirs that are in 
contact with the qubits, I express the pumping and decay rates 
in terms of new parameters Ti and n ll49ll : 



ri = ri{i-ri), Pi = riri (/ = i,2). 



(3) 



The range < r^ < 1, includes a medium that only absorbs 
excitation (decay, r^ = 0), one which only provides it (pump, 
r/ = 1), as well as a the most common assumption of a ther- 
mal bath with finite temperature, (white noise, r^ < 1/2). The 
parameter Ti = Yi-\- Pi quantifies the interaction of each qubit 
with its reservoir as well as its effective spectral broadening. 

Taking into account the evolution of the density matrix p of 
our bipartite system, one can conclude that in the steady state, 
it has the general block diagonal form: 



/Poo 

Pii P12 

Pi*2 P22 

V P33 



with 



Poo = 1 - (^1) - (^2) + (^1^2) , 
Pii = {ni)-{n\n2), /= 1,2, 

Pl2 = (^12)*, 
P33 = (^1^2) , 



(4) 



(5a) 
(5b) 
(5c) 
(5d) 



in terms of the operators rii = G- Gi and n\2 = gIg2. The av- 
erage value {fii) is the probability for qubit / to be excited. 



regardless of the other one. (^12) accounts for the popula- 
tion transfer between the qubits and (^1^2) for their effec- 
tive coupling, in the sense that, independent qubits would lead 
to (^1^2) = (^i)(^2). The general expressions for the steady 
state of two coupled qubits can be written as {nt) = Pf^ /Tf^, 

(ni«2) = (/'i{«2> + /'2(«i))/(ri + r2), (ni2) = '-^M^^ 

where the effective parameters Pf^ = Pi + {P\ -\- P2)Xi and 
r?^ = Fi + (Fi + r2)Xi are expressed in terms of an effec- 
tive coherent exchange factor Xt 



V/rs,-,! 



the Purcell rates 14711 . Xi quantifies how efficiently the exter- 
nal inputs and outputs are distributed among the qubits thanks 
to the coherent coupling and despite the total decoherence, 

Hot = ri +12 + rf + T^. 

In Ref. ll47n . I showed that the steady state p of two cou- 
pled qubits is the same than that of a four-level system, that 
is, the system depicted in Fig.[Hb) with no correspondence to 
two two-level systems in Fig. [Ha), but rather to a single en- 
tity. This is the case of four single atomic levels or of a single 
quantum dot that can host two excitons and form a biexciton 
state. The results presented in the following sections are di- 
rectly based on p or on averaged (single-time) quantities (O) 
computed as Tr(pO) and, therefore, are also valid for a four- 
level system. 



III. ENTANGLEMENT AND LINEAR ENTROPY 

Among the four Bell states, \(j)±) and | V^±), only the last 
two are achievable in the present configuration (since po3 = 0). 
Let us therefore write the most general entangled state that can 
be achieved in this system as | i/a) = (| 1) +^'^ | 2))/V^. Log- 
ically, the larger the probability to find the qubits in such state 
(the closer p is to | i/a) (i/A |), the larger is the degree of entan- 
glement in the mixture represented by p. In order to make 
this statement mathematically precise, we can make explicit 
the entangled contribution to p by expressing it as 

P=Poo|0)(0|+p33|3)(3| 

+ 7?i|1)(1|+7?2|2)(2|+7?^|Va)(va| (6) 

where Ri = pn - |pi2|, R2 = P22 - IP12I andT?^, = 2|pi2|. Rt 
(/ = 1 , 2, i/A) are not probabilities (Ri,R2 may be negative) but, 
when they are normalized as 



Ri 



\Ri 



P00 + P33 + |^l| + |^2|+^V.' 



(7) 



they represent the contribution of the pure states | /) to the 
mixture where the entangled state has been identified and set 
apart. In order to enhance entanglement, we must maximize 
Ryr (pn) while minimizing the populations poo and P33, and 
the differences ^1 and ^2- The non-entangled contributions 
can be put together in a single expression to be minimized: 

^=1-V 

The degree of entanglement can be quantified by the con- 
currence (C) (501], which ranges from (separable states) 
to 1 (maximally entangled states). It defined as C = 



[max{0, vAi — v^ — v^ — v^}] , where { Ai , A2, A3 , A4} 
are the eigenvalues in decreasing order of the matrix 
pTp*T, with T being a anti diagonal matrix with elements 
{ — 1,1,1,-1}. The concurrence in this system is given by 
C = 2Max[{0, IP12I — ^/poop33}]^ which shows a threshold 
behaviour that we anticipated above: the coherence between 
the intermediate states, |pi2|, must overcome the population 
of the spurious states poo, P33. A related important factor 
to build some concurrence, is the degree of purity in the 
system [i5 lil . This is measured through the linear entropy, 
Sl = ^[1 — Tr(p^)], which is for a pure state, and 1 for a 
maximally disordered state (where all four states occur with 
the same probability 1/4). 

Without loss of generality, we can analyze the entanglement 
and linear entropy of our system in the steady state, by consid- 
ering the parameters A > 0, on the one hand, < r2 < ri < 1 
on the other hand, and arbitrary Fi , F2 > 0. This simply im- 
plies that we label as 2 the qubit that is in contact with the 
medium which has the most dissipative nature. Let us ig- 
nore dephasing effects for the moment (we bring them back 
in Sec. El). 

We start by noting that if ri = r2 = r, that is, if the reser- 
voirs are of the same nature, there is no entanglement in the 
system (C = 0), regardless of all the other parameters, since, 
in this case, the density matrix is diagonal with elements 
{(1 — r)^,r(l — r),r(l — r),r^}, that is, a mixture of separable 
states. This result has already been pointed out in the litera- 
ture 1251 127|], however, let us insist on the fact that, it is not 
the amount of decoherence induced on the qubits by their en- 
vironments what destroys entanglement, but their similarity in 
nature (or temperatures in the case of thermal baths). Let us 
then consider the cases with r2 < ri in the rest of this section. 

As in Ref. Dili , I examine the region of the concurrence- 
linear entropy plane that our system can access in Fig. [SJa). 
The shaded region is reconstructed by randomly choosing all 
the parameters and computing their C and Sl. The accessi- 
ble region is well below the black thin line for the maximally 
entangled mixed states Dili , that provides the maximum con- 
currence achievable for a given linear entropy. More inter- 
estingly, the points are bounded in good approximation by a 
second (dashed blue) line specific to our system. This line 
corresponds to the extreme case of reservoirs with exactly op- 
posite natures, ri = 1 and r2 = 0, but equally strong influence 
on the qubits, Fi = F2 = F (that is, Pi = 72 = F, ^2 = 71= 0). 
The steady state can be written in terms of a single unit-less 
complex number, 

ae'^ = {^-iT)/g, (8) 

with norm a and phase j3: poo = P22 = P33 = 4:^, Pn = 



4+a2 



and P12 



4+a2 



The two qubits are sharing a single 



excitation (^^ ) + (^2) = 1 • Note that both detuning (A) and the 
average decoherence (F) contribute symmetrically to a and 
have the same effect on the steady state: to make the coherent 
coupling less effective. The phase j3 = — arctan (F/A), which 
is the same than that of the entangled state | i/a) formed in the 
steady state, can be rotated by changing these two parameters. 
This is a way to phase shift the entangled state obtained in the 
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system. Concurrence and linear entropy read 
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FIG. 2: (a) Distribution in the C-Si plane of all the possible two 
qubit configurations (shaded region). The thin solid line corresponds 
to the maximum C for a given Si in a general bipartite system. The 
dashed blue line corresponds to the optimal configuration (ri = \, 
r2 = 0, Fi = r2 = F, see Eq. (l9at), a good approximation to the 
maximal C vs Si in our system (with the exception of the dark blue 
region above). Below, in dark purple, the particular case of thermal 
baths, (b) C (solid bl ack) and Si (dashed red) for the optimal case as 
a function of a = VA^ + F^/g. In inset (c), the non-entangled con- 
tributions to the steady state: ^i (dotted brown), R2 (dashed purple) 
and R (solid blue). The shaded area represents Ry/ (as R^-^-R = 1). 
Idem in (d) and (e), but for two thermal baths at infinite and zero 
temperatures: ri = 1/2, r2 = 0. The vertical guide lines mark the 
points where entanglement appears and where it is maximum, close 
to the points where R is minimum and R approaches ^1 , respectively. 
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}], (9a) 



(9b) 



I plot them in Fig. [2b), as a function of a moving leftwards 
along the dashed blue line of Fig. [2a). The contributions ^1 
(dotted brown), R2 (dashed purple), R (solid blue) and ^1^ 
(shaded area) are also presented in the inset. Fig. [2c), for a 
better understanding of the origin of entanglement. At van- 
ishing a (or large effective coupling), the excitation is equally 
shared among the states and the system is maximally mixed 
(with Ry^ = 0). Concurrence becomes different from zero at 
a = 1, which is close to the point where the non-entangled 
contribution to the density matrix, R reaches its minimum (and 
Ry/ its maximum). The contributions of the spurious states 
Poo, p33 have been considerably reduced while the coherence 
IP12I is sufficiently large to overcome them. 
The maximum concurrence in the absolute for this system. 



Cu 



5-l)/4^31%, 



(10) 



is reached at a = 1 + >/5. This is the region where a large 
contribution Ry/ is combined with low Si. Moreover, the non- 
entangled contribution R becomes similar to ^1 meaning that 
the steady state is close to the mixture Myf = T?,^ | i/a)(i/a | 
-\-{l—Ry/) I 1)(1 |, only between the entangled state | Xj/) and 
I 1). The large contribution of | 1) to the steady state is ex- 
pected since the first qubit is pumped and the second decays. 
What is less expected is that, by populating this state, we are 
purifying the total mixture and enhancing the presence of the 
entangled state. Increasing a further leads to the saturation 
of the system into state My/ and eventually to self-quenching 
of coherence D2I] . Note, however, that concurrence decreases 
slowly and never becomes strictly zero again, due to the fact 
that p -^ Myf and, therefore, C -^ Ry/. 

The small region in dark blue in Fig. [2a), to the right and 
above the line in Eq. ([9aj) , corresponds to cases more entan- 
gled for the same entropy, than the configuration previously 
discussed. Relaxing the previous conditions to Fi 7^ r2, for 
instance, is enough to fill this area. In any case, configura- 
tions above the dashed line exist only for very mixed states, 
with Sl > (17 — 3\/5)/30 ^ 0.34, being less appealing for 
applications. In Fig. [3a), I plot the corresponding concur- 
rence as a function of Fi and F2 in order to show that it is 
robust to their difference: C > as long as Fi + F2 > 2, and 
C > Cniax/2 ~ 15% in most of the area shown. 

To conclude this analysis, one can check that if one medium 
provides an overall dissipation and the other one, an overall 
gain (0 < r2 < 1/2 < ri < 1), then C can reach non-negligible 
values (above 10%). This is one of the important results in this 
text, the opposite nature of the reservoirs can lead the steady 




FIG. 3: Contour plots of concurrence C as a function of Fi and F2 
for (a) ri = 1, r2 = and (b) ri = 1/2, r2 = 0. To the left of the 
white lines, C = 0. The maximum value achieved with this system is 
in (a), Cmax ~ 31%. In (b), with thermal baths, concurrence rises up 
to C ^ 10% for an asymmetric configuration. 



state close to M,^, allowing for the highest degrees of entan- 
glement in the system. 

On the other hand, if we keep equal interactions with the 
reservoirs, Fi = r2 = F, but they are not restricted in their 
natures, the values of the concurrence decreases. Let us con- 
sider the case that has been previously studied in the literature, 
two thermal reservoirs in contact with two qubits. In our nota- 
tion: < r2 < ri < 1/2. In this case, the concurrence does not 
grow higher than 4% (as we said in the introduction), which is 
reached for the extreme case, ri = 1/2, r2 = 0. Again, within 
the new restrictions, it is favorable for entanglement that ex- 
citation is provided through one qubit while the other only 
dissipates. The equivalent expressions to Eqs. Q read: 



C = Max[(«,^^^Mp!Zl„, 



Sl 



4 + a2 
3(4 + a2)2 • 



(11a) 
(lib) 



This case is featured in Fig.[2d)-(e). We observe that C be- 
comes different from zero again at the minimal R (maximal 
entangled contribution Ryf) and that its maximum value is 
reached when ^1 approaches R. However, the concurrence 
remains one order of magnitude smaller than Cmax and the 
linear entropy does not drop to 0. With thermal excitation, Ry/ 
is always too small and the steady state is not close enough 
to My/ to exhibit a high degree of entanglement. However, 
in contrast with the optimally pumped case, one can increase 
entanglement from these figures by allowing Fi 7^ F2. Con- 
currence is increased, filling the purple darker shaded region 
in Fig. [3a). The maximum concurrence here is C ^ 10% at 
Fi ^ 1.24 and F2 ^ 6.45. This is shown in Fig. Ob) where the 
highest values of C appear in light grey around those rates. 



IV. ANTIBUNCHING 

Is there an experimental observable that can evidence the 
high degrees of entanglement that we have analysed here? 



One possibility is to reconstruct the steady state density ma- 
trix through quantum tomography, but this method involves 
complicated set ups and numerous and repeated measure- 
ments [53]. Here, I propose an alternative method that only 
involves photon counting, that is, the quantity: 

5 = {ni){n2)-{nin2) = Pnp22- Poop33' (12) 

(^1) and (^2) are proportional to the intensity of the light 
emitted from each qubit, obtained by counting photons from 
each source, while (^1^2) is obtained by counting simulta- 
neous photon detections. 5 is directly linked to the second 
order cross correlation function ll47ll at zero delay, ^^2 (^) ~ 
1 — 5/{{n\){n2))- S is zero if the qubits, acting as two ran- 

^2) (2) 

dom photon sources, are independent (g\2 (0) = g\2 (°q) = 1), 
and different from zero if one qubit 's emission is conditional 
to the other qubit's state. 5 < implies that the simulta- 
neous emission from both qubits is enhanced in the system 
as compared to the independent emissions. This is a neces- 
sary condition for photon bunching, although bunching also 

requires ^^2 (^) -^ ^12 v'^)- ^^ example where 5 < 0, is the 
cross simultaneous emission of two coupled harmonic oscilla- 
tors 11461] . On the other hand, 5 > implies that simultaneous 
emission from both qubits is less likely than in the uncoupled 
situation. Again, this is necessary for photon antibunching 

In the steady state of our system, the emission from one 

(2) (2) 

of the qubits is always antibunched (g-- ^(0) = < g-- \t), 
/ = 1,2, as it correspond to a two-level system) and the cross 
emission from both qubits fits < 5 < 1/4. One can check 
these limits from gathering 5 s from many randomly gener- 
ated configurations. It cannot go below zero because the 
only coherence and entanglement in the system come from 
the state | Xj/) (that gives the maximum value 5 = 1/4) and 
not I 0) = (I 0) + e'P' I 3))/V^ (that would give 5 = -1/4). 
The sign of 5 is linked to the type of entangled state realised in 
the system, also when there is superposition or mixture with 
other states and C < 1 . For instance, if weplot C versus 5 
for the maximally entangled mixed states Dili with entangle- 
ment provided by | i/a), we obtain the black thin line in Fig.JH 
that abruptly falls at 5 = 1/9. This is because for this kind of 
states with C < 2/3, 5 remains constant. One would obtain 
the symmetrical curve at negative 5, if the mixture was with 
I 0). In this example, high degrees of entanglement are related 
to large 5. The other dotted lines appearing in Fig. [Hare more 
examples of this relationship between the type of entangle- 
ment and C with 5. The central black dotted line corresponds 
to the superposition or a mixture of | Xj/) with | 0), when | Xj/) 
is the dominant state: C = 45. There is a symmetric counter- 
part curve (not shown) in the opposite situation, where | 0) is 
dominant, with C = 4| 5 1 and 5 < 0. The upper red dotted line 
corresponds to the superposition or mixture of | Xj/) with | 0) 
or with I 3): C = 2v5. The counterpart curve, with | 0), is 
symmetrical. The space between these two dotted lines could 
be filled with mixtures of | Xj/) with both | 0) and | 3). The 
lower blue dotted line corresponds to the mixture of | Xf/) with 
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1 — Vl —45. In all these cases, large 
5 is correlated with large C, although the connexion is rather 



the state My^: C = 
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FIG. 4: Distribution in the C-d plane of all the possible qubit con- 
figurations (shaded region). The thin solid and dotted lines corre- 
spond to different examples of entangled mixed states, with | Xj/) 
(see the main text). The pure state | \j/) corresponds to the extreme 
point (1/4, 1). The dashed blue line corresponds to the configuration 
(ri = l,r2 = 0, Fi = r2 = F, increasing F anti-clockwise) which en- 
closes all possible realizations. Inside, in dark purple, the particular 
case of thermal reservoirs. In inset: C vs (^i) for thermal reservoirs 
(where (ni) < 0.5, in dark purple) and for those configurations where 
5 > 0.04 (in blue). The dashed blue line (ri = l,r2 = 0, Fi = F2 = F) 
goes clockwise with increasing F. 



general and not exclusive enough to define any entanglement 
witness in terms of 5. 

Let us go back to our system and investigate how to use 
these correlations to extract information about C from the 
measured 5. The shaded region in Fig. |4l corresponds, as in 
Fig. [21 to the situations realised in our system. It is completely 
enclosed this time by the dashed blue line, which corresponds 
to reservoirs with opposite natures (as analysed in the previous 
section). In this limiting case, 5 reads 



V4 + a2y 



(13) 



Thanks to this analytical boundary, we can turn the general 
statement that there is some correlation between 5 and C into 
a more accurate (mathematical) one: C-(5) < C < C+(5) 
where 



C±(5) 



2V25\/l ± Vl - 165 - 85 - V28 
l±Vl-165 



(14) 



These inequalities become most stringent when 5 > 0.04, for 
instance 5 > 0.061 implies 20% < C < 28.3%. More precise 
information can be obtained if (^1) is included in the analy- 
sis, looking at the inset of Fig. IH In blue, we see a cloud of 
numerically generated points where 5 > 0.04. There is also 
a clear correlation between large, unsaturated, population of 
the dot (0.8<(^i)<0.91) and large degrees of entanglement 

(10%<C<Cn,ax). 

Such large 5 and C, cannot be obtained with thermal reser- 
voirs for the qubits. The small accessible area in that case, 
is shaded in darker purple in Fig. |4] (within the dashed blue 
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FIG. 5: Effect of dephasing on the results for opposite reservoirs 
{ri = \, r2= 0, Fi = F2 = F), of C (a) and 3 (b) as a function of 
F. The set of curves corresponds to values of y^ from to 20g, 
increasing in steps of 2g. Entanglement (a) is diminished by pure 
dephasing, going from top to bottom curves, and the maximum is 
reached at higher Fs (see the dashed line, joining the maxima of all 
curves). The maximum 3 remains 1/16 for all values of dephasing 
although this is reached for lower Fs. 



curve) and in the inset. One cannot use 5 or (^1) as entangle- 
ment indicators because, for all 5 or (^1), C can take a broad 
range of values that always includes 0. 



V. PURE DEPHASING 

Pure dephasing provides extra decay for the coherence in 
the system. It weakens the correlations established between 
the qubits and is, therefore, an enemy of entanglement. We 
can see this in Fig. Oa) where I plot the effect of increasing 
dephasing {yf = y^ = y^) on entanglement. The curve at the 
top is the same as in Fig. [2b), with opposite kinds of reser- 
voirs. The rest of the curves correspond to increasing values 
of the dephasing rate in steps of 2g up to y^ = 20g. Entan- 
glement decreases, and its maximum value for a given y^, re- 
quires higher F. The set of maximum C and the corresponding 
required F are plotted with a dashed thin line superimposed to 
the curves for clarity. Entanglement is quite robust in this con- 
figuration, it disappears but asymptotically and very slowly. 
Note that y^ must be one order of magnitude larger than g so 
that C is decreased to the values obtained with thermal reser- 
voirs (10%). 

In Fig. Ob), I plot the counterpart curves for 5, that shrink 
and move leftwards with dephasing. However, the maximum 
5 remains 1/16 for all dephasing, taking place at lower Fs. 
Given that the tendency of the maximum 5 is the opposite 
to that of the maximum C, the possibility of using 5 as an 
indicator of entanglement fails at large dephasing. 



VI. CONCLUSIONS 

I have computed the entanglement (C) and linear entropy 
{Sl) for two coupled qubits in the steady state created by an 
incoherent continuous excitation. I have studied, not only the 
case where the excitation is of a thermal origin, but also a 
more general out-of-equilibrium situation where populations 
can be inverted ((^i) > 0.5). In this case, I find that entan- 
glement can be greatly enhanced (C up to 31%) as compared 
to the best thermal values (C up to 10%), with also a much 



higher purity of the state. This is obtained in a configura- 
tion where one qubit essentially dissipates excitation while 
the other essentially gains it. I have used both numerical re- 
sults, in the most general case, and analytical formula, in this 
optimal case, to fully understand and characterise entangle- 
ment formation. Entanglement (provided by the entangled 
state I V^)) is enhanced in the steady state when the pumped 
qubit approaches saturation, because this removes population 



from spurious states (| 0) and | 3)) and purifies the statistical 
mixture, even in the presence of pure dephasing. Finally, I 
have shown that the quantity 5 = (^i)(n2) — (^1^2), that can 
be measured experimentally by photon counting, can be used 
as an indicator of high degrees of entanglement in this system 
and specially for the optimal configuration. 

Acknowledgments: This research was supported by the 
Newton International Fellowship program. 



[1 
[2: 
[3 

[4: 

[5 
[6 

[7 
[8 
[9 

[lo: 

[11 
[12: 

[13 

[14: 

[15 

[16: 

[17 
[18 

[19: 

[20 
[21 

[22 

[23 
[24: 

[25 
[26: 

[27 

[28 

[29: 



R. Horodecki, P. Horodecki, M. Horodecki, and K. Horodecki, 

Rev. Mod. Phys. 81, 865 (2009). 

M. A. Nielsen and I. L. Chuang, Quantum computation and 

quantum information (Cambridge University Press, 2000). 

T. D. Ladd, F. Jelezko, R. Laflamme, Y. Nakamura, C. Monroe, 

and J. L. OBrien, Nature 464, 45 (2010). 

G. W. Gardiner and P. Zoller, Quantum Noise (Springer- Verlag, 

Berlin, 2000), 2nd ed. 

Z. Ficek and R. Tanas, Phys. Rep. 372, 369 (2002). 

F. Verstraete, M. M. Wolf, and J. I. Cirac, Nature Phys. 5, 633 

(2009). 

T. Yu and J. H. Eberly, Science 323, 598 (2009). 

Z. Ficek, Front. Phys. China 5, 26 (2010). 

D. Braun, Phys. Rev. Lett. 89, 277901 (2002). 

M. S. Kim, J. Lee, D. Ahn, and R L. Knight, Phys. Rev. A 65, 
040101(R) (2002). 

L. Jakbczyk, J. Phys. A: Math. Theor. 35, 6383 (2002). 
S. Schneider and G. J. Milburn, Phys. Rev. A 65, 042107 
(2002). 

F. Benatti, R. Floreanini, and M. Piani, Phys. Rev. Lett. 91, 
070402 (2003). 

L. Xiang-Ping, F. Mao-Fa, Z. Xiao- Juan, and C. Jian-Wu, Chi- 
nese Phys. Lett. 23, 3138 (2006). 

E. del Valle, F. P. Laussy, and C. Tejedor, Europhys. Lett. 80, 
57001 (2007). 

E. del Valle, F. P. Laussy, F. Troiani, and C. Tejedor, Phys. Rev. 
B 76, 235317 (2007). 

L. D. Contreras-PuHdo and R. Aguado, Phys. Rev. B 77, 
155420 (2008). 

M. Hor-Meyll, A. Auyuanet, C. V. S. Borges, A. Aragao, 
J. A. O. Huguenin, A. Z. Khoury, and L. Davidovich, Phys. 
Rev. A 80, 042327 (2009). 

F. Benatti, R. Floreanini, and U. Marzolino, Phys. Rev. A 81, 
012105 (2010). 

L. Jakbczyk, J. Phys. B: At. Mol. Opt. Phys. 43, 015502 (2010). 

X. X. Yi, C. S. Yu, L. Zhou, and H. S. Song, Phys. Rev. A 68, 

052304 (2003). 

M. B. Plenio and S. F Huelga, Phys. Rev. Lett. 88, 197901 

(2002). 

X. Hao and S. Zhub, Eur. Phys. J. D 41, 199 (2007). 

X. L. Huang, J. L. Guo, and X. X. Yi, Phys. Rev. A 80, 054301 

(2009). 

J.-B. Xu and S.-B. Li, New J. Phys. 7, 72 (2005). 

L. Hartmann, W. Diir, and H.-J. Briegel, Phys. Rev. A 74, 

052304 (2006). 

L. Hartmann, W. Dur, and H. J. Briegel, New J. Phys. 9, 230 

(2007). 

N. Lambert, R. Aguado, and T. Brandes, Phys. Rev. B 75, 

045340 (2007). 

A. Rivas, N. P. Oxtoby, and S. F. Huelga, Eur. Phys. J. B 69, 51 



(2009). 
[30] H. Wang, S. Liu, and J. He, Phys. Rev. E 79, 041113 (2009). 
[31] L. Zhou, G. H. Yang, and A. K. Patnaik, Phys. Rev. A 79, 

062102 (2009). 
[32] J. Li and G. S. Paraoanu, New J. Phys. 41, 113020 (2009). 
[33] C.-J. Shan, T. Chen, J.-B. Liu, W.-W. Cheng, T.-K. Liu, Y.-X. 

Huang, and H. Li, Int. J. Theor. Phys. 49, 717 (2010). 
[34] I. Bloch, Nature 453, 1016 (2008). 
[35] M. Bayer, P. Hawrylak, K. Hinzer, S. Fafard, M. Korkusinski, 

Z. R. Wasilewski, O. Stern, and A. Forchel, Science 291, 451 

(2001). 
[36] H. J. Krenner, M. Sabathil, E. C. Clark, A. Kress, D. Schuh, 

M. Bichler, G. Abstreiter, and J. J. Finley, Phys. Rev. Lett. 94, 

057402 (2005). 
[37] B. D. Gerardot, S. Strauf, M. J. A. de Dood, A. M. Bychkov, 

A. Badolato, K. Hennessy, E. L. Hu, D. Bouwmeester, and P. M. 

Petroff, Phys. Rev. Lett. 95, 137403 (2005). 
[38] J. Clarke and F. K. Wilhelm, Nature 453, 1031 (2008). 
[39] A. Imamoglu, D. D. Awschalom, G. Burkard, D. P. DiVincenzo, 

D. Loss, M. Sherwin, and A. Small, Phys. Rev. Lett. 83, 4204 

(1999). 
[40] S. B. Zheng and G. C. Guo, Phys. Rev. Lett. 85, 2392 (2000). 
[41] S. Ashhab, A. O. Niskanen, K. Harrabi, Y Nakamura, T. Picot, 

P C. de Groot, C. J. P M. Harmans, J. E. Mooij, and F Nori, 

Phys. Rev. B 77, 014510 (2008). 
[42] S. Osnaghi, P. Bertet, A. Auffeves, P. Maioh, M. Brune, J. M. 

Raimond, and S. Haroche, Phys. Rev. Lett. 87, 037902 (2001). 
[43] J. Majer, J. M. Chow, J. M. Gambetta, J. Koch, B. R. Johnson, 

J. A. Schreier, L. Frunzio, D. I. Schuster, A. A. Houck, A. Wall- 
raff, et al.. Nature 449, 443 (2007). 
[44] A. Laucht, N. Hauke, J. M. Villas-Boas, F Hofbauer, 

M. Kaniber, G. Bohm, and J. J. Finley arXiv:0904.4759 (2010). 
[45] E. Gallardo, L. J. Martinez, A. K. Nowak, D. Sarkar, H. P. 

van der Meulen, J. M. Calleja, C. Tejedor, I. Prieto, D. Grana- 

dos, A. G. Taboada, et al., Phys. Rev. B 81, 193301 (2010). 
[46] E. del Valle, Microcavity Quantum Electrodynamics (VDM 

Verlag, Saarbriicken, Germany, 2010). 
[47] E. del Valle, Phys. Rev. A 81, 053811 (2010). 
[48] H. J. Carmichael, Statistical methods in quantum optics 1 

(Springer, 2002), 2nd ed. 
[49] H.-J. Briegel and B.-G. Englert, Phys. Rev. A 47, 3311 (1993). 
[50] W K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 
[51] W J. Munro, D. F. V. James, A. G. White, and R G. Kwiat, 

Phys. Rev. A 64, 030302(R) (2001). 
[52] O. Benson and Y Yamamoto, Phys. Rev. A 59, 4756 (1999). 
[53] M. Steffen, M. Ansmann, R. C. Bialczak, N. Katz, E. Lucero, 

R. McDermott, M. Neeley, E. M. Weig, A. N. Cleland, and J. M. 

Martinis, Science 313, 1423 (2006). 



